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Turbulent  Transport  at  High 
Reynolds  Numbers  in  an  Inertial 
Confinement  Fusion  Context 

Mix  is  a  critical  input  to  hydro  simulations  used  in  modeling  chemical  or  nuclear  reac 
tion  processes  in  fluids.  It  has  been  identified  as  a  possible  cause  of  performance  degra 
dation  in  inertial  confinement  fusion  (ICF)  targets.  Mix  contributes  to  numerical  solution 
uncertainty  through  its  dependence  on  turbulent  transport  coefficients,  themselves  uncer 
tain  and  even  controversial  quantities.  These  coefficients  are  a  central  object  of  study  in 
this  paper,  carried  out  in  an  Richtmyer  Meshkov  unstable  circular  two  dimensional  (2D) 
geometiy  suggested  by  an  ICF  design.  We  study  a  pre  turbulent  regime  and  a  fully  devel 
oped  regime.  The  former,  at  times  between  the  first  shock  passage  and  reshock,  is  charac 
terized  by  mixing  in  the  form  of  interpenetrating  but  coherent  fingers  and  the  latter,  at 
times  after  reshock,  has  fully  developed  turbulent  structures.  This  paper  focuses  on  the 
scaling  of  spatial  averages  of  turbulence  coefficients  under  mesh  refinement  and  under  var 
iation  of  molecular  viscosity  [i.e.,  Reynolds  number  (Re)].  We  find  that  the  coefficients 
scale  under  mesh  refinement  with  a  power  of  spatial  grid  spacing  derived  from  the  Kolmo 
gorov  213  law,  especially  after  reshock.  We  document  the  dominance  of  turbulent  over  mo 
lecular  transport  and  convergence  of  the  turbulent  transport  coefficients  in  the  infinite  Re 
limit.  The  transport  coefficients  do  not  coincide  for  the  pre  and  post  reshock  flow  regimes, 
with  significantly  stronger  transport  coefficients  after  reshock.  [DOI:  10.1 1 15/1.4027382] 

Keywords:  computational  fluid  dynamics,  hydrodynamics,  large  eddy  simulation  (LES), 
mixing,  turbulence  modeling,  shocks 


1  Introduction 

Turbulence,  starting  with  Re  ~  2000  [1],  profoundly  alters  the 
fluid  state.  Such  flows  show  a  large  increase  in  fluid  transport, 
known  as  turbulent  transport.  Mix  (either  molecular  level  mix 
ing,  or  fine  scale  granularity  of  flow  mixtures)  can  occur  in  the 
absence  of  turbulence.  The  rate  of  mixing  however,  is  greatly 
enhanced  by  the  presence  of  turbulence.  Mix  of  reactive  compo 
nents  can  enhance  reactions,  while  admixture  of  nonreactive  con 
stituents  can  retard  it.  Thus,  we  see  the  importance  of  turbulent 
mix  to  theories  of  turbulent  combustion  and  turbulent  nuclear  re 
active  processes.  We  analyze  flows  having  Re  fa  3.5  x  104  and 
higher;  the  starting  value  is  set  somewhat  arbitrarily  from  the 
Reynolds  numbers  achieved  in  laboratory  experiments,  cf.  [2],  of 
the  Rayleigh  Taylor  (RT)  instability.  The  appeal  to  RT  experi 
ments  is  based  on  the  strong  sensitivity  of  the  RT  macroscopic 
observables,  i.e.,  the  overall  mixing  rates,  to  turbulent  transport 
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[3,4].  In  addition,  when  analysis  consists  of  quantities  subject 
to  numerical  uncertainties  and  even  controversies,  comparison 
to  experiment  (validation)  is  an  indispensable  strategy.  A  simu 
Iation  validation  protocol  based  on  these  ideas  is  described 
in  [5], 

The  purpose  of  this  paper  is  to  explore  the  dependence  of  turbu 
lent  transport  parameters  on  mesh  and  Re,  for  a  range  of  Re 
including  ICF  applications  and  a  range  of  mesh  which  includes 
the  finest  three  dimensional  (3D)  meshes  currently  feasible  for  a 
complex  engineering  problem  of  this  nature.  This  study,  with  its 
multiple  simulations  and  variation  of  parameters,  is  not  practical 
to  carry  out  in  3D  and  so  we  regard  this  two  dimensional  (2D) 
study  as  offering  insight  to  possible  future  3D  simulations.  In  con 
trast  to  3D  turbulence,  2D  turbulence  has  two  conserved  quantities 
(energy  and  enstrophy).  The  turbulent  energy  spectrum  has  two 
asymptotic  scaling  ranges,  a  classical  Kolmogorov  scaling  k~5/3 
for  intermediate  wave  numbers  k  and  a  specifically  2D  scaling 
law  k~3  driven  by  enstrophy  for  large  k  up  to  the  dissipation 
range.  This  theory  was  presented  in  the  classical  arguments  of 
Kraichnan  [6]  and  observed  in  direct  numerical  simulations 
(DNS)  of  Boffetta  [7].  The  energy  flows  in  opposite  directions 
under  these  two  cascades,  with  their  associated  flow  length  scales, 
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while  the  wave  number  for  energy  injection  marks  the  division 
between  them. 

The  context  selected  for  this  study  is  a  Richtmyer  Meshkov 
(RM)  fluid  instability  in  a  circular  geometry.  Aside  from  its  inter 
est  as  a  classical  hydrodynamic  instability,  this  problem  is  sugges 
tive  of  scientific  issues  which  arise  in  studies  of  ICF  capsules. 
Three  stages  in  an  ICF  implosion  are  subject  to  mix  and  resulting 
performance  degradation.  These  are  the  (i)  ablation  stage  (RT 
instability  modified  by  mass  transfer  effects),  (ii)  shock  passage 
through  various  capsule  layers  (RM  instabilities),  and  (iii)  decel 
eration  stage,  which  leads  to  hot  spot  formation  (RT  instability). 
The  role  of  mix  in  each  of  these  stages  is  the  subject  of  ongoing 
studies.  The  present  study  pertains  to  the  second  stage,  shock  pas 
sage,  with  a  possible  relevance  to  the  other  ICF  mix  stages. 

We  study  two  qualitatively  different  flow  regimes,  the  pre  and 
post  reshock  flows.  The  pre  reshock  regime  is  incipiently  turbu 
lent  and  characterized  by  mixing  at  the  level  of  inteipenetrating 
but  coherent  fingers.  The  post  reshock  regime  is  characterized  by 
fully  developed  turbulence  and  a  proliferation  of  fine  scale 
mixing. 

In  a  high  Re  shock  driven  flow,  molecules  ionize  and  exhibit 
strong  thermal  coupling  to  the  electron  and  radiation  fields.  For 
this  reason,  we  consider  a  fluid  with  a  plasmalike  highly  heat  con 
ductive  Prandtl  number  Pr  =  10~4  and  a  Schmidt  number  Sc  =  1. 
These  are  laminar  transport  values.  The  primary  purpose  of  this 
paper  is  to  examine  the  corresponding  turbulent  transport  coeffi 
cients.  We  define  the  high  Re  limit  as  molecular  viscosity  v  — >  0 
with  fixed  values  for  Sc  and  Pr.  This  theoretical  definition  yields  a 
uniquely  specified  theoretical  limit.  To  complete  the  specification 
of  fluid  transport  input  parameters,  we  set  the  (molecular  level) 
isotropic  viscosity  to  zero  and  introduce  a  Re  dependent  shear 
(anisotropic)  kinematic  viscosity  i/  =  5.13x  10-4  m2/s  for 
Re  =  6  X  105.  We  use  inverse  scaling  of  v  with  Re  for  other 
Reynolds  numbers.  The  levels  of  mixing  and  the  value  of  Re 
achieved  in  ICF  processes  is  a  subject  of  ongoing  research,  out 
side  of  the  scope  of  this  study.  We  study  a  series  of  mesh  levels 
(Ay  from  6.25  x  10~4  to  1.56  x  lCT4m)  and  Re  values  from 
Re  =  3.5  x  104  to  a  finite  mesh  approximation  of  Re  ~  oo 
(obtained  by  setting  the  molecular  viscosity  to  0).  We  mention 
previous  studies  of  the  same  problem  [5,8].  The  numerical  study 
is  based  on  the  front  tracking  code  FronTier  [9,10]  enhanced  with 
an  LES  turbulence  model  (Germano  and  Moin’s  dynamic  subgrid 
model  [11,12])  and  dynamic  subgrid  scale  (SGS)  terms.  The 
tracking  reduces  excess  numerical  species  diffusion. 

The  SGS  terms  are  defined  using  filters  and  test  filters.  For  each 
term  we  introduce  a  model  or  approximate  form  with  a  to  be 
determined  coefficient.  The  filters  are  convolution  operators 
which  average  the  numerical  solution  with  approximate  localiza 
tion.  We  choose  a  unit  mesh  block  average  for  the  filter.  The  test 
filter  is  one  unit  coarser  in  resolution;  for  this  we  choose  a  2 D 
array  of  mesh  blocks  in  D  dimensional  space.  The  basic  idea 
expressed  in  Germano’s  identity  [11]  is  that  the  differences 
between  the  SGS  stresses  at  the  test  filter  level  and  the  filter  level 
can  be  evaluated  (due  to  cancellation)  in  terms  of  quantities 
defined  at  a  mesh  level  and  otherwise  missing  coefficients  can  be 
determined.  This  computable  difference  of  stresses,  known  as  the 
Leonard  term  [13],  is  denoted  as  L. 

Suppressing  some  details  [14],  the  SGS  model  cM  is  a  vector 
quantity  whose  divergence  appears  in  the  RHS  of  the  conservation 
law.  It  has  the  form  cM  =  cAv2p|S|V{7,  where  Ay  is  the  mesh 
spacing,  S  is  the  strain  rate  matrix,  and  VU  is  a  solution  gradient. 
[In  the  various  conservation  laws,  VU  becomes  the  anisotropic 
(isotropic)  part  of  the  strain  rate  matrix  for  the  anisotropic  (iso 
tropic)  turbulent  flux,  the  temperature  gradient  multiplied  by  a 
specific  heat  for  the  thermal  flux  or  the  species  concentration  gra 
dient  for  species  flux.]  The  dimensionless  coefficient  c  is  to  be 
determined  separately  for  each  of  the  conservation  equations;  the 
rest  of  this  expression  is  the  model  M. 

Let  Mdlff  denote  the  difference  between  the  filter  dependent 
models  Mflller  and  M,estfil,er  for  the  SGS  terms  taken  at  the  filter 
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and  test  filter  grid  levels.  Setting  the  difference  in  models 
(assumed  to  have  a  common  coefficient)  proportional  to  the  differ 
ences  of  stresses  L,  we  have  the  equation 

L  =  cMdiff  (1) 

where  the  same  coefficient  c  relates  models  to  unclosed  SGS 
terms  at  the  filter  and  test  filter  levels.  The  key  equation  is  then 
c  =  L/Mdlff.  The  kinematic  turbulent  transport  coefficient  j;turb  is 
given  by 

Xturb  =  cAy2^!  (2) 

with  c  dimensionless  and  S  the  strain  rate  matrix.  Here  Xturb 
denotes  one  of  the  turbulent  shear  viscosity,  turbulent  bulk  viscos 
ity,  turbulent  species  diffusion,  or  turbulent  thermal  diffusion,  in 
each  case  with  a  distinct  value  for  c.  The  dimensions  of  j;turb  are 
[l2/t]. 

The  form  of  the  SGS  model  M  =  fourb pVU  as  a  product  of  two 
first  order  spatial  derivatives  multiplied  by  Ay2  suggests  that  spa 
tial  averages  of  xturb  should  be  governed  by  Kolmogorov’s  2/3 
law,  scaling  as  Ay4/3.  We  see  this  relation  definitively  after 
reshock  and  a  weak  association  at  best  to  this  relation  before 
reshock.  The  scaling  holds  not  only  for  the  velocity  field  but  for 
the  concentration  and  thermal  fields  as  well. 

In  Sec.  2,  we  examine  the  turbulent  energy  spectrum  for  our 
simulations.  In  Sec.  3,  we  relate  the  scaling  laws  for  spatial  mean 
values  (defined  by  a  simple  average  over  the  mixing  zone)  of  the 
turbulent  transport  coefficients  to  scaling  laws  for  Kolmogorov 
spectra.  Using  the  scaling  laws,  we  study  the  mean  values  of  the 
transport  coefficients  as  Ay  — >  0  and  Re  ~  oo.  Contrasts  are 
drawn  between  the  singly  shocked  regime  (before  reshock),  which 
is  pre  turbulent  and  the  doubly  shocked,  post  reshock  regime  with 
fully  developed  turbulence.  Conclusions  are  discussed  in  Sec.  4. 


2  Scaling  Laws  for  Primitive  Equation  Quantities 

To  examine  the  scaling  properties  of  the  turbulent  spectral 
energy,  we  plot  in  Fig.  1  log  log  graphs  of  the  turbulent  energy 
spectrum  E(k)  versus  wave  number  k  for  the  Re  ~  oo  data  taken 
from  the  finest  simulation  performed.  The  spectral  analysis 
depends  on  a  Fourier  expansion.  In  view  of  the  annulus  shape  of 
the  mixing  zone,  we  analyze  with  the  approximate  bounding  box 
for  the  mixing  zone.  This  results  in  a  wave  number  normalized  by 
size  of  the  bounding  box  around  the  mixing  zone  (approximately 
3/8  of  the  domain  size).  The  left  frame  shows  data  from  before 
reshock  and  the  right  frame  after  reshock.  In  the  left  frame  we  see 
a  small  wave  number  spectral  range  for  A~5/3  and  beyond  this,  a 
steeper  decay  less  than  the  theoretical  value  of  3  but  possibly 
influenced  by  this  effect.  The  turbulence  in  the  left  frame  is  not 
yet  fully  developed  at  this  time  and  thus  the  flow  features  are  out 
side  of  the  scope  of  Kraichnan’s  theory.  This  imperfect  scaling 
can  also  be  observed  in  Table  1.  In  the  right  frame  (after  reshock), 
Fig.  1  suggests  two  scaling  regimes,  an  approximate  A~5/3  and  an 
approximate  A'-3,  which  is  consistent  with  Kraichnan’s  theory. 
This  change  in  scaling  law  can  also  be  observed  in  Table  2  with 
grid  levels  I  II  showing  better  agreement  with  each  other  than 
with  grid  level  III.  The  Kolmogorov  scale  for  the  flows  we  com 
pute  varies  from  microns  or  less  to  tens  of  microns  at  the  time  of 
reshock,  while  the  finest  mesh  considered  is  Ay  =  1.56  x  10~4  m. 
The  possible  A~3  scaling  regime  is  visible  for  only  1/2  an  order  of 
magnitude  in  wave  numbers  and  thus  is  not  definitive,  which  is  in 
line  with  comments  from  [7],  Beyond  the  differences  in  physical 
modeling  between  [7]  and  the  present  work,  we  do  not  come  close 
to  resolving  the  Kolmogorov  scale.  For  a  compressible  calcula 
tion,  we  are  not  remotely  near  to  DNS.  Our  finest  mesh  spacing 
Ay  =  1.56  x  10~4  m  is  about  a  factor  of  100  larger  than  our  esti 
mated  Kolmogorov  scale  ij  ~  L  ■  Re~3/4  ~  2.5  x  10  6m,  based 
on  the  length  L  =  5.5  x  10~2  m  (width  of  the  mixing  zone)  and 
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Fig.  1  Log-log  energy  spectrum  plots  versus  wave  number.  The  reference  lines  have  slopes  k5!3  (dashed)  and  k  3  (dashed- 
dotted).  Left:  Pre-reshock.  Right:  Post-reshock. 


Table  1  Scaled  mean  turbulent  transport  coefficients.  Mesh  I:  400x800,  with  mesh  II  doubled  and  mesh  III  doubled  again.  Units 
are  m2/s  before  scaling.  Time  t  5.8x10  5s,  midway  between  the  first  and  second  shock  passage.  In  the  column  Re  ss  the  value 
is  obtained  by  setting  the  molecular  viscosity  to  0. 


Mesh 

Re  3.5  x  104 

Re  6  x  105 

Re  6  x  106 

Re  6  x  107 

Re  «  oo 

cS'  t^rb / Ax4/3  scaled  anisotropic  turbulent  viscosity  (m2'3/s) 

i 

23 

24 

24 

24 

24 

ii 

19 

20 

20 

20 

21 

hi 

14 

14 

14 

16 

15 

C^turb/ At4' 3  scaled  isotropic  turbulent  viscosity  (m2'3  / s) 

i 

86 

102 

95 

77 

87 

ii 

65 

81 

75 

65 

64 

hi 

51 

62 

60 

54 

50 

cS'  fi turb/Av4/3  scaled  species  turbulent  diffusion  (m2/3/s) 

i 

37 

52 

49 

38 

42 

ii 

29 

39 

37 

29 

28 

hi 

21 

26 

28 

23 

20 

cS'  atUrb / Ax4/3  scaled  turbulent  thermal  diffusion  (m2/3  / s) 

i 

103 

71 

44 

40 

47 

ii 

65 

66 

42 

33 

34 

m 

30 

56 

41 

33 

28 

Table  2  Scaled  mean  turbulent  transport  coefficients.  Meshes  and  units  as  in  Table  1.  Time  f  =  9.0  x  10  5  s,  shortly  after  the  pas¬ 
sage  of  reshock,  in  the  column  Re  «  the  value  is  obtained  by  setting  the  molecular  viscosity  to  0. 


Mesh 

Re  3.5  x  104 

Re  6  x  105  Re  6  x  106 

Re  6  x  107 

Re  ~  oo 

I 

44 

cS1 

vL,/A*4/3 

51 

scaled  anisotropic  turbulent  viscosity  (m2/3/s) 

53 

47 

54 

II 

48 

48 

56 

58 

53 

III 

57 

49 

64 

64 

62 

I 

214 

cS1 

Vturb/A*4/3 

248 

scaled  isotropic  turbulent  viscosity  (m2/3/s) 

238 

214 

234 

II 

222 

236 

234 

232 

235 

III 

253 

233 

276 

277 

260 

I 

112 

cS ' 

A urb/A*4/3 

136 

scaled  species  turbulent  diffusion  (m2/3/s) 

127 

127 

124 

II 

115 

133 

128 

122 

124 

III 

133 

122 

149 

148 

133 

I 

126 

cS’ 

aturb/A.v4/3 

163 

scaled  turbulent  thermal  diffusion  (m2/3/s) 

143 

122 

132 

II 

106 

181 

158 

135 

145 

III 

103 

174 

191 

174 

149 
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Fig.  2  Fractional  mesh  error  for  dimensionless  turbulent  transport  coefficients,  comparing  coarse  to  fine  (l-lll)  and  medium  to 
fine  (ll-lll)  grids,  plotted  versus  Re.  Left:  Before  reshock.  Right:  After  reshock.  Curves  labeled  1,  2,  3  denote  inverse  isotropic 
viscosity,  Schmidt  number,  and  Prandtl  number,  respectively.  The  dashed-dotted  line  denotes  the  error  in  the  comparison  l-lll 
and  the  solid  line  denotes  the  error  in  the  comparison  ll-lll. 


Fig.  3  Re  dependence  of  mean  turbulent  transport  coefficients  for  a  Richtmyer-Meshkov  instability.  Fractional  variation  for 
each  of  the  four-dimensional  transport  coefficients,  plotted  as  xturb(Re)//turb(Re  sa  «)  versus  Re,  using  the  finest  grid  level. 
Left:  Before  reshock.  Right:  After  reshock.  Curves  labeled  0,  1,  2,  3  denote  anisotropic  viscosity,  isotropic  viscosity,  species 
diffusivity,  and  thermal  diffusivity  shown  as  a  fraction  of  the  values  of  these  parameters  at  Re  «  »  and  plotted  versus  the 
Reynolds  number. 


Re  =  600, 000.  This  simple  analysis  suggests  many  unexplored 
length  scales  between  the  mesh  resolution  and  the  Kolmogorov 
scale,  a  statement  which  remains  true  even  for  the  smallest 
Re(Re  =  35,000)  considered  here. 

The  statistical  nature  of  the  space  time  dependent  pdfs  varies  as 
the  distinct  scaling  ranges  (intermediate,  large  k,  dissipation 
range)  change.  Thus,  mesh  convergence  requires  that  the  coarse 
graining  length  scale  [15  18]  used  to  define  the  pdfs  must  be  cho 
sen  to  lie  within  a  single  scaling  range  (if  the  mesh  is  chosen  finer, 
i.e.,  not  in  this  range,  as  with  Ax  — >  0). 

3  Convergence  and  Variation  of  Mean  Turbulent 
Transport  Coefficients 

Assuming  that  c  is  mesh  convergent,  the  turbulent  transport 
coefficient  ^turb  defined  by  (2)  cannot  be  asymptotically  independ 
ent  of  Ay.  The  strain  rate  matrix  S  is  assumed  to  satisfy  the 
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Kolmogorov  2/3  law  for  the  local  behavior  of  the  two  point  corre 
lations,  leading  to  a  bound  proportional  to  Ay  ' 2' 3  for  each  deriva 
tive  and  for  each  turbulent  transport  term,  and  A*-4' 3  for  S  factors 
combined  with  the  gradient,  part  of  M  in  the  flux  term.  We  rewrite 
(2)  to  obtain 

*turb  =  cAy4/3  (W/3)  =  cAx^S'  (3) 

where  (3)  defines  S'.  Since  c  is  assumed  to  converge  to  a  finite 
limit  as  Ay  — >  0  and  S'  is  assumed  to  be  bounded  due  to  a  Kolmo 
gorov  type  scaling  law,  we  see  that  /tarb  vanishes  with  Ay  at  a  pre 
dieted  rate.  Remarkably,  we  observe  the  same  decay  rate  for  all  of 
the  turbulent  coefficients,  even  though  the  decay  law  is  derived 
from  an  incompressible,  velocity  correlation  decay  (the  Kolmo 
gorov  2/3  law  for  the  decay  of  the  velocity  correlations).  This 
scaling  law  is  supported  for  the  fully  developed  turbulence  of 
Table  2  and  is  marginal  at  best  for  Table  1.  It  is  clear  that  these 
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Reynolds  Number 


Fig.  4  Turbulent  transport  as  a  fraction  of  total  transport  plotted  versus  Re  for  each  of  four  mean  transport  coefficients,  with 
data  taken  at  the  fine  grid  level.  Left:  Before  reshock.  Right:  After  reshock.  Curves  labeled  0, 1 , 2,  3  denote  anisotropic  viscosity, 
isotropic  viscosity,  species  diffusivity,  and  thermal  diffusivity,  respectively,  as  a  fraction  of  total  transport. 


two  flow  regimes,  even  in  the  infinite  Re  limit,  define  qualitatively 
different  turbulent  flows.  Additionally,  the  post  reshock  coefficient 
mean  values  show  only  mild  finite  Re  effects  for  Re  >  3.5  x  104, 
while  the  pre  reshock  coefficients  transition  to  (different)  nearly  Re 
independent  values  at  higher  values  of  Re.  The  turbulent  thermal 
and  species  coefficients  are  nearly  identical  to  each  other. 

In  Fig.  2,  we  plot  fractional  mesh  errors  versus  Re,  comparing 
coarse  to  fine  and  medium  to  fine,  for  the  three  dimensionless 
mean  turbulent  transport  coefficients.  The  formula  used  for  calcu 
lation  of  this  error  is  | f  /|/|/'|,  where/  stands  for  the  mean  turbu 
lent  dimensionless  coefficient  computed  at  the  finest  grid  and  / 
stands  for  the  corresponding  parameter  at  a  medium/coarse  grid 
level.  We  see  that  the  mesh  errors  are  either  small  in  relative  terms 
or  convergent  under  mesh  refinement  or  both  for  all  values  of  Re. 
There  is  an  exception  for  the  turbulent  Prandtl  number  and 
high  Re  turbulent  Schmidt  number,  especially  before  reshock, 
for  which  further  mesh  refinement  studies  would  be  helpful. 
Further  levels  of  mesh  refinement  would  allow  a  deeper  exam 
ination  of  2D  turbulence  scaling  properties  for  the  present 
problem. 

To  better  understand  the  high  Re  asymptotics,  we  plot  in  Fig.  3 
the  fractional  change  ;(turb(Re)/;(turb(Re  ~  °°)  in  the  mean  dimen 
sional  transport  coefficients  versus  Re  at  the  finest  grid.  The  dif 
ferences  between  the  pre  and  post  reshock  turbulent  flow  can  be 
seen  by  comparing  the  left  frame  to  the  right.  We  observe  a  rather 
mild  Re  dependence,  especially  for  the  fully  developed  turbulent 
flow  (right). 

To  see  the  dominance  of  the  turbulent  transport  as  a  fraction  of 
the  total  transport,  we  plot  Xturb/Xtotai  versus  Re  in  Fig.  4.  Aside 
from  the  Prandtl  number,  the  Re  ~  oo  values  are  nearly  reached 
by  Re  =  6  x  106  for  the  pre  turbulent  flow  (left)  and  by 
Re  =  6  x  105  for  fully  turbulent  flow  (right).  Note  the  striking 
qualitative  difference  between  the  Prandtl  number  and  the  other 
transport  coefficients.  Due  to  the  very  large  level  of  laminar  (mo 
lecular)  heat  conductivity  assumed,  even  the  high  but  finite 
Re  =  6  x  107  simulations  display  strong  finite  Re  number  effects 
for  Prtotai.  It  is  for  this  reason  that  the  Prandtl  number  curves  do 
not  follow  the  pattern  set  by  the  other  transport  coefficients.  The 
convergence  of  the  total  transport  to  its  Re  ~  oo  asymptotic  val 
ues  is  slower  for  the  pre  turbulent  (singly  shocked)  flow  than  it  is 
for  the  fully  developed  turbulent  (doubly  shocked)  flow,  a  conse 
quence  of  the  smaller  values  of  turbulent  transport  for  the 
pre  turbulent  flow  (singly  shocked). 
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4  Conclusion 

We  have  shown  the  mesh  dependence  for  the  spatially  averaged 
turbulent  transport  coefficients  to  be  either  relatively  small  or  con 
vergent  or  both,  for  all  values  of  Re  considered  and  for  the  range 
of  length  scales  (Ax)  considered.  The  notion  of  pdf  convergence 
depends  on  an  observational  length  scale  for  coarse  graining  simu 
lation  values  to  define  a  local  pdf.  Convergence  is  limited  to  a 
range  of  length  scales  for  which  constant  scaling  behavior  is  valid. 
The  post  reshock  coefficients  have  a  mild  Re  dependence,  while 
finite  Re  effects  in  these  coefficients  persist  to  higher  Re  values 
for  the  pre  turbulent,  single  shocked  case.  The  common  post 
reshock  Kolmogorov  2/3  law  scaling  among  these  coefficients  is 
unexplained.  The  high  level  of  thermal  conductivity  leads  to  finite 
Re  effects  in  the  thermal  transport  even  for  quite  elevated  values 
of  Re. 
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